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Exact solutions with an exponential behaviour of the scale factors are considered in a multidimensional cosmological 
model describing the dynamics of n + 1 Ricci-flat factor spaces Mi in the presence of a one-component perfect 
fluid. The pressures in all spaces are proportional to the density: pi — wtp, i = 0, n. Solutions with accelerated 
expansion of our 3-space Mo and a small enough variation of the gravitational constant G are found. These solutions 
exist for two branches of the parameter wo . The first branch describes superstiff matter with wo > 1 , the second 
one may contain phantom matter with wo < — 1. 



1. Introduction 

Fundamental physical constants, relations between them 
and their possible variations are a reflection of the situ- 
ation with unification [HEHHIUEIE] . 

Here we are mainly interested in the gravitational 
constant and its possible variations. The oldest prob- 
lem is that of possible temporal variation of G, which 
arose due to papers by Milne (1935) and Dirac (1937). 
In Russia, these ideas were developed in the 60s and 
70s by K.P. Staniukovich [7][5], who was the first to 
consider simultaneous variations of several fundamental 
constants. 

Our first calculations based on general relativity with 
a perfect fluid and a conformal scalar field [8] gave G/G 
at the level of 10 -11 — 10~ 13 per year. Our calcula- 
tions in string- like [9] and multidimensional models with 
a perfect fluid [10] gave the level 10~ 12 , those based 
on a general class of scalar-tensor theories [13] and a 
simple multidimensional model with p-branes [121114] 
gave for the present values of cosmological parameters 
10~ 13 — 10~ 14 and 10~ 13 per year, respectively. Similar 
estimations were made by Miyazaki within Machian the- 
ories [15] giving for G/G the estimate 10~ 13 per year 
and by Fujii — on the level 10~ 14 — 10~ 15 per year [IB] . 
Analysis of one more multidimensional model with two 
curvatures in different factor spaces gave an estimate on 
the level 10 -12 [17]. Here we continue our studies of 
variation of G in one more multidimensional cosmolog- 
ical model with perfect fluid. 

2. The model 

We consider a cosmological model describing the dy- 
namics of n Ricci-flat spaces in the presence of a one- 
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component "perfect-fluid" matter 
the model 



The metric of 



g = - exp[2~f(t)]dt <g> dt + exppa^i)]^ 



i=0 



is defined on the manifold 
M = Rx M a x ...x M n , 



(1) 



(2) 



where Mi with the metric g % is a Ricci-flat space of di- 
mension di, i = 0, . . . , n; n > 2. The multidimensional 
Hilbert-Einstein equations have the form 

1 



— 8ff R — K Tjy , 



(3) 



where k 2 is the gravitational constant, and the energy- 
momentum tensor is adopted as 



(T N ) = diag(-p, p S™°, p n S^ n n ), 



(4) 



describing, in general, an anisotropic fluid. 

We assume the pressures of this "perfect" fluid in all 
spaces to be proportional to the density, 



Pi (t) = (l-Ui/di) P (t), 



(5) 



where m = const , i = 0, . . . , n . We also put p ^ . 

We also impose the following restriction on the vec- 
tor u = (ui) € W n+1 : 



(u,«),^0. 

Here, the bilinear form 
defined by the relation 

(u, v )* = G lJ UiVj, 

u, v 6 R™ +1 , where 

8 i: > 1 

G %3 = - — h 

<&• 2-D 



tn+l 



pn+1 



(6) 
I is 

(7) 
(8) 



are components of the matrix inverse to the matrix of 
the minisuperspace metric [2!?ll2"3"] 



Gij — diSij — didj. 



(9) 
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In (|S]), D = 1 + 5Z^L dj is the total dimension of the 
manifold M ©. 

The restriction ([5]) reads 



i=0 



i=0 



3. Solutions with exponential scale 
factors 

Here, we consider a special family of solutions with 
an exponential behaviour of the scale factors from [T8l 
119) with the metric written in the synchronous time 
parametrization 



-dt s ® dt s + ^2af(t s )g\ 



(11) 



i=0 



Solutions with a an exponential behaviour of the scale 
factors take place for 



(u w -u,u)* = 0. 
Here and below the vector 



v^=2d i 



(12) 



(13) 



corresponds to the A-term "fluid" with pi = —p (vacuum- 
like matter). 

In this case, the solutions are determined by the met- 
ric (TTTj) with the scale factors 



a L = a,i(t s ) = A 4 exp(z/i s ), 
and the density 

p = const. 
Here 



eu a — 



2k 2 P 



(14) 
(15) 

(16) 



where e = ±1, u % = G^Ui and Ai are positive con- 
stants, i = 0,...,n. Here p > for (it, u)* < and 
p < for (u,u)* > 0. 

The model under consideration was integrated in [18] 
for (u, it)* < 0. The solutions from [TH] were gener- 
alized in 19J to the case when a massless minimally 
coupled scalar field was added. Families of exceptional 
solutions with power-law and exponential behaviours of 
the scale factors in terms of synchronous time were sin- 
gled out in pjj] and correspond to a constant value of 
the scalar field: (p — const. When the scalar field is 
omitted, we are led to solutions presented in [4] and 
above for power-law and exponential cases, respectively 
(in [18] these solutions were originally written in the har- 
monic time parametrization) . It may be verified that the 
exceptional solutions with an exponential dependence 
of the scale factors are also valid when the restriction 
< is replaced by ©. 



4. Acceleration and variation of G 

In this section, the metric g° is assumed to be flat, 
and d = 3. The subspace (M ,g°) describes "our" 
3-dimensional space and (Mi^g 1 ) the internal factor 
spaces. 

We are interested in solutions with accelerated ex- 
pansion of our space and small enough variations of the 
gravitational constant obeying the present experimental 
constraints, see 1121: 



|G/(GiT)|(* s o)<0.1, 
where 

«o 



H = 



(17) 



(18) 



is the Hubble parameter. We suppose that the internal 
spaces are compact. Hence our 4-dimensional constant 
is (see [TU]) 



G = const ■ n"_ 1 ( a i~ d * 



(19) 



We will use the following explicit formulae for the 
contravariant components: 



3=0 

and the scalar product reads 
(u^ -u,u), 



.1 1 



(20) 




(21) 



4.1. Exponential expansion with acceleration 

For solutions under consideration, an accelerated expan- 
sion of our space takes place for 



v° >0. 



(22) 



Remark: D= 4 case. For D = 4, when the internal 
spaces are absent, we get u° = —ug/6 and 



,( A ) 



1 



u, w), = - (u - 6)u = 0, 
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(23) 
(24) 



which implies Uq = 6, or, equivalently, 

p = -p. (25) 
We get 



(26) 



which agrees with the well-known result for D — A: the 
de-Sitter solution with the cosmological constant A = 
K 2 p > . The condition h>° > is equivalent to e = — 1 . 
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For our exponential solutions we get 



3=1 

and hence 



n 

G/(GH) = --J2^di 



= 5. 



(27) 



(28) 



3=1 



The constant parameter S describes variation of the 
gravitational constant and, according to (fT7j) ■ 

|<5|<0.1. (29) 

It follows from the definition of v % in (|T6|) that 



1 " 



or, in terms of covariant components (see l|20p). 



6= - 



(D-4)uo-2YZ =1 Ui 
1(5- D)u + Y%=i u i' 



(30) 



(31) 



Thus the relations (J2TJ), dHj) , J29J), (|3T|) and the con- 
straint (fT2|) determine a set of parameters Ui compatible 
with the acceleration and tests on G-dot. 

In what follows we will show that these relations do 
really determine a non-empty set of parameters Ui de- 
scribing the equations of state. 

4.1.1. The case of constant G 

Consider the most important case (5 = 0, i.e., when the 
variation of G is absent: G = 0. 

Indeed, there is a tendency of lowering the upper 
bound on G . Moreover, according to arguments of [21], 
S < 10~ 4 . This severe constraint just follows from the 
identity 



G/G = a/a 



(32) 



that takes place in some multidimensional models. Here 
a is the fine structure constant. 

Isotropic case. We first consider the isotropic case 
when the pressures coincide in all internal spaces. This 
takes place when 

ui = vdi, i = 1, ... ^n. (33) 

For pressures in internal spaces we get from ([5]) 

Pi = (l-v)p, i=l,...,n. (34) 

Then we get from ((lOj) and (j2T|) 



,.( A ) 



1 



2 - D 



{-Ud-l)u + 2du v-2dv 2 }, (35) 



u, u)* = 



1 



■ [2u + 2dv 



2- D 1 

+ h(d- l)ul - 2du v + 2dv 2 ] 



(36) 



Here and in what follows we denote d = D — 4 . 
For 6 = , we get in the isotropic case 

v = uo/2, (37) 

or, in terms of pressures, 

Pi = (3p - p)/2, 2 = 1,..., n. (38) 

Substituting (|37|) into (|35|) and (f36|) . we get 

(u,u)* = -ul/6, (39) 
(u (A) -u,u)* =ii (u -6)/6. (40) 

Here, we obtain the same relations as for D = 4. 
For our solution, we should put ito 7^ and hence, due 
to d]), 

(41) 



tt = 6, 
i.e., 
Po = —pi 



Pi 



-2p, 



i > 0. 



(42) 



Using ([33]) and ([37]) . we get u° = — uq/6 = —1 and 
= for i > 0, hence z/j = for i = 1, . . . , n, i.e., all 
internal spaces are static. 

The metric (fTTj) reads in our case 

n 

g = -dt s ® cft s + exp(2^ t s )<7 + ^ ^.g 1 , (43) 



;.=i 



where Ai are positive constants, and 



2 



(44) 



For accelerated expansion we get e = — 1 . We see that 
the power v° is the same as in case D = 4 . 

Anisotropic case. Consider the anisotropic (w.r.t. in- 
ternal spaces) case with 5 = 0, or, equivalently (sec 

eu), 



(45) 



(D - A)u = 2^2ui. 

i=l 

This implies 

( m (A) - it, u )* = ^uo("o - 6) - A, 
( it, u )* 
where 



±7/ 2 

6 U 



(46) 
(47) 



" 2 1 / n 



i — 1 x z — 1 



> 0, d = D — A. (48) 



The inequality in (|48[) can be readily proved using 
the well-known Cauchy-Schwarz inequality: 



> 



(49) 



Indeed, substituting = \fd~i and Cj = Ui/^fdi into 
(gSJ, we get (gSJ). The equality in (gl]) takes place only 



4 



when the vectors (bi) and (cj) are linearly dependent, 
which for our choice reads: Ui/^fdi = vy/dl where v is 
constant. Thus A = only in the isotropic case 
In the anisotropic case we get A > . 

In what follows we will use the relation 



(u (A) - u,u)* = i(it - u£)(u - u ), 
where 



it 







3± V9 + 6A 



(50) 
(51) 



are roots of the quadratic trinomial (I46p obeying 

Uq < 0, u+ > 6 for A > 0. (52) 

-ito/6 and hence 



It follows from (1451) that u° 
u 



12n 2 p 



6 V u 2 , - 6A ' 



Here (see ((T 



U = Uq. 



(53) 



(54) 



Accelerated expansion of our space takes place when 
v a > 0, or, equivalently, when either 



(A) uq = u , e = 1 or 

(B) U = Uq , e = -l. 

In terms of the parameter wq . 
Po = wop, w = 1 - u /3, 
these two branches read: 
(A) w = w- = Jl+ |A, 



(B) Wo = w+ = - A /l 



A. 



(55) 
(56) 



(57) 

(58) 
(59) 



The first branch (A) describes super-stiff matter ( wo > 
1 ) with negative density p < . 

The second branch (B) corresponds to matter with 
positive density (since (u, u)* < 0). This matter is 
phantom (i.e., wo < —1) when A > 0. 

4.1.2. The case of varying G 

Now we consider the case S ^ 0, i.e., when G / 0. 
In what follows, we use the observational bound (f29|) : 
| S | < 0.1, stating the smallness of S. 
Using (|3"Tj) . we get 



i 



-dbuo, 



where d = D — 4 and 

l + <5(l-c/)/(3d) 



6 = 6(5) = 



1 — 5/2 



(60) 



(61) 



For the scalar product we get from (fBT)]) 



u, m )* = i Au 2 , — i?u — A, 



,(A) _„ „\. = 1 
6 ' 

>,u)* = -^Au 2 , + A ' 



where A was defined in (J48J) (see <JT0]) and (f2~T]) ). 

A 1 / d,\ 2 d, 2 1 
6 rf + 2\ 2/4 3 



B = 



1 



d + 2 



(2 + cffi). 



Using (|6ip . we obtain the explicit formulae 
(d + 2)<5 2 



A = A(S) = 1 
5 = B{8) = 



12d(l- J/2) 2 ' 
1 -5/3 



1 -(5/2' 



(62) 
(63) 



(64) 
(65) 

(66) 
(67) 



Due to |(5| < 0.1, A is positive, A > 0, and close to 
unity: |A - 1| < ±10" 2 . 

For the contravariant component u° we get from 
(1201) and (il: 



(68) 
(69) 



w° = -Cuo/6, 
where 

C = C(<5) = 3B - 2 = 1/(1 - 5/2). 
It follows from (63) and J68|) that (see ([16])) 



n Cwo / 12k 2 p 

v = — e 



6 V Au 2 - 6A ' 



Here (due to (fT2"|) ) 



u = u 



±(5) = -AZB ± ^9S 2 + 6AA) 



A 



(70) 



(71) 



are roots of the quadratic trinomial (f6"2")l . 

Isotropic case. Let us consider the isotropic case (|33]). 
Then we obtain from (l60|) 

v = dbuo/2, (72) 

or, in terms of pressures 

Pi = ^[3bp + (2-3b)p], i=l,...,n. (73) 



For scalar products we get 

= -Auq/6, 
(m (A) -u,u)* = u (Aii - 65)/6. 



(74) 
(75) 



For our solution, we should put «o ^ and hence 
Mo = 6B/A > 0. The metric (ITT]) reads in our case 

n 

g = -dt. ® dt s + Al e 2 ^g° + e 2 ^ £ A 2 g\ (76) 
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where A4 are positive constants, 



Cu 12k 2 p 



and 



(77) 



V = V = £ 



Su 



12k 2 p 
Any 



1,. 



U{l-5/2) 

n . The last formula follows from (| 16[) and 
u S 



(78) 



6d(l - 5/2)' 



(79) 



We see that the power v° does not coincide, for <5 7^ 0, 
with that in case D = 4 . 

The accelerated expansion condition for our 3D 
space, v° > , reads in this case 



u 



6.8(5) 



-1 



A(S) ' 

or, equivalently, in terms of w ([57} (po = ^op) 
2B(5) 



(80) 



wo = w£ (8) = 1 



A(S) 



For <5 > , we get isotropic contraction of the whole 
internal space Mi x . . . x M n . In this case, Wq (8) < — 1 , 
and hence phantom matter occurs with the equation of 
state close to the vacuum one since 



w+(8) + l 



5(1 + S/d) 



3[l-5 + (d-l)S 2 /(6d)Y 



(82) 



For small 6 we have w£(6) = -1 - 8/3 + 0(S 2 ) . 

For 8 < we obtain isotropic expansion of the whole 
internal space. Then Wq (8) > —1 , and phantom matter 
does not occur. 

Anisotropic case. Now we consider the anisotropic 
case A > when 6^0. Here uq — u (S), sec (fTTj) . 

Accelerated expansion of our 3-dimensional space 
takes place when v > 0, or, equivalently, when either 



(A) u = u (8), 

(B) u = u+(5), 



£ = 1 
£ = -1. 



or 



(83 
(84 



1 



In terms of the parameter wq (po 
- mq/3 ) these two branches read 



w p, w 



(A) w 

(B) w 

where 



w a (8), 

+ 




w+(8), 



wt(8) = l-4(8)/3. 
For small 8 we have 

w b k (S)=w±(0)-~(l± 



(85 
(86 

(87 



V9 + 6A 



0{8 2 ). 



Thus for small 8 the parameter wq has a small de- 
viation from that obtained for 8 = 0. For small 8 , wq 
shifts by 0(8) term. 



The first branch (A) describes superstiff matter, 
wq > 1 , since (8) > 1 due to Uq (8) < . It may 
be shown that the density is negative in this case since 
(u,u)* > 0. 

For branch (B) we find that 



w+(S) < -1 
only if 

A > 6[A{8) - B{8)} 



-8/{l-8/2f. 



(89) 



(90) 



This is a condition on the appearance of "phantom" 
matter. For 8 > this inequality is valid, but for 
8 < it is satisfied only for a big enough value of the 
anisotropy parameter A , see (|90| . 

5. Conclusions 

We have considered multidimensional cosmological mod- 
els describing the dynamics of n + 1 Ricci-flat fac- 
tor spaces Mi in the presence of a one-component 
anisotropic (perfect) fluid with pressures in all spaces 
proportional to the density: pi = Wip, i = 0, ...,n. So- 
lutions with accelerated expansion of our 3-dimensional 
space Mq and small enough variation of the gravita- 
tional constant G were found. In the general non- 
isotropic case, these solutions exist for two branches of 
the parameter wo: (A) from ([55]) and (B) from (f8"6"]) . 
Branch (A) describes super-stiff matter with wo > 1 
while branch (B) may corresponds to phantom matter 
with wo < — 1 • 
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